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ABSTRACT 

The partial widths are determined for squark decays to gluinos and quarks, and gluino 
decays to squarks and quarks, respectively. The widths are calculated including one-loop 
SUSY-QCD corrections. The corrections amount to +30% to +50% for squark decays and 
— 10% to +10% for gluino decays. We have derived the results in the DR and MS renor- 
malization schemes, and we have demonstrated explicitly that the one-loop effective qqg and 
qqg couplings are equal in the limit of exact supersymmetry. 
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1 Introduction 



The theoretical predictions for the decay properties of squarks and gluinos in supersymmetric 
theories are of interest for several reasons. The decay properties determine the experimental 
search techniques, and the measurement of decay branching ratios will allow us to study 
the couplings between the standard and the novel particles, interrelated by supersymmetry. 
A thorough analysis of the partial widths including higher-order corrections is required for 
this purpose at the theoretical level. SUSY-QCD corrections, which take into account the 
spectrum of the standard quarks/gluons and their supersymmetric partners squarks/gluinos, 
have so far been performed only for a few examples: The production of squark and gluino 
pairs at hadron colliders |]] and of squark pairs in e + e~ annihilation 0, the decay of squarks 
to quarks and photinos 0, and the decay of charged Higgs bosons to stop/sbottom || and 
top/bottom particles In this note we present the theoretical predictions for the decay 
channels 

q — ► g + q for nig > rrig (1) 

g — > q + q / q + q for m~ g > rriq (2) 

including SUSY-QCD corrections. The analysis is performed for light quarks, excluding the 
case of top (s) quarks. This calculation is more demanding than the calculation for the decay 
processes discussed earlier in the literature due to the three colored particles involved in the 
initial and final states of @j and @). At the technical level, we will follow the path traced 
out for the production of squarks and gluinos in hadron collisions [|TJ. 

2 Theoretical Set-up 

In analyzing the SUSY-QCD corrections, the supersymmetric squarks are assumed to be 
degenerate in mass. The formulae for squark decays apply to both chiral states, qi and (Jr. 
For gluino decays they include the sum over all flavor and chiral components and the two 
possible final states in Eq. (Q) related by C conjugation. The final states involving top quarks 
are not taken into account since we have treated the quarks in the final states as massless 
particles; the case of (s)top particles will be deferred to a subsequent analysis. In loops we 
have used a top-quark mass of m t = 180 GeV and we have taken the top-squarks to be 
mass-degenerate with the other squarks. The masses of the squarks and gluinos are assumed 
to be above 100 GeV. (For experimental bounds see Ref. 0). The squark and gluino masses 
are defined as pole masses. 

The diagrams for the squark and gluino decays are displayed in Fig. |T[ While (a) repre- 
sents the Born diagrams for the two decay modes, the self-energy diagrams are shown in (b), 
vertex corrections in (c), and hard-gluon radiation in (d). In (c) and (d) we only give the 
squark-decay diagrams, as the gluino and squark decay processes are related to each other 
by crossing. 
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The ultra-violet, infrared, and collinear divergences have been regularized by switching 
to n = 4 + e dimensions. The actual calculation has been carried out in two schemes: 

(i) It is well-known that dimensional reduction DR M preserves supersymmetry at least 
to one-loop. This scheme is hence an obvious candidate for performing the renormalization; 
some technical difficulties in handling the factorization of collinear singularities, however, 
still demand a proper understanding ||. 

(ii) From a technical point of view, the MS renormalization scheme || is easier to handle 
in general. However, by introducing a mismatch between the number of gauge-boson and 
gaugino degrees of freedom in n ^ 4 dimensions, the scheme violates supersymmetry explic- 
itly. In particular, the qqg Yukawa coupling g s , which by supersymmetry should coincide 
with the qqg gauge coupling g s , deviates from g s by a finite amount at the one- loop order. 
Requiring the physical amplitudes to be independent of the renormalization scheme, a shift 



between the bare Yukawa and gauge couplings must be introduced in the MS scheme [10 



9s = 9s 



1 + 



(3) 



which effectively subtracts the contributions of the false, non-supersymmetric degrees of 
freedom (also called e scalars); Ca = 3 and Cf = 4/3 are the Casimir invariants of the 
SU (3) gauge group. 



Taking into account this shift in the MS scheme, we find that all results are identical 
with the results obtained within the DR scheme bearing in mind the relation between the 
gauge couplings in the two schemes. This comparison provides non-trivial consistency checks 
of the one-loop SUSY-QCD corrections. 

The need for introducing a finite shift is best demonstrated for the effective qqg Yukawa 
coupling, which must be equal to the qqg gauge coupling in a supersymmetric world with 
massless gauginos/gluons and equal- mass squarks/quarks. For the sake of simplicity we 
define the effective couplings g e J^{Q 2 ) and g e s ^{Q 2 ) in the limit of on-shell squarks/quarks 
and almost on-shell gluinos/gluons with virtuality Q 2 <C m? = m 2 . [In this limit the 
couplings do not involve gauge- dependent terms.] 

First we consider the corrections of the couplings 

g e s ff (Q 2 )=g s (l + S) and gff(Q 2 )=g s (l + 5) (4) 



in the MS scheme. By adding the contributions of the external self-energies and the charge 
renormalization to the vertex corrections [all given explicitly in the appendix], we find the 
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following results[|: 



MS : 5(qqg) 



5(qqg) 



2 1 / Q 2 
^2 l0g Ul 



+ 



C, 



OCa 

Air 



2 i /g 2 

e 2 \m~ 



+ 



6 



(5) 
(6) 



The difference between 5 and 5 is indeed compensated by the finite shift introduced in 
Eq. (|). As a result, the couplings g e J^ (Q 2 ) and g e J^ (Q 2 ) agree with each other to one-loop, 
as required by supersymmetry. The same observation can be made in the DR scheme. In 
this scheme we find 



DR 



8{qqg) 



5{qqg) 



47T 



T^ Ca 



2 1 1 



2 1 i 
log 

£ 2 & 



I 



mi 



(7) 



As anticipated, the DR renormalization scheme preserves supersymmetry and guarantees 
that the gauge and Yukawa couplings also coincide to one-loop order in the supersymmetric 
limit. (Note the difference between the effective gauge couplings in the MS and DR scheme, 
(a s /4vr) x C A /6; see also Ref. [TO, IT1). 



3 Results 

At the Born level, the partial widths of the decay processes ([!]) and (0) are given by 

Uq^~9 + q) = ^ { Q 3 9 (9) 
2 m~ 

a e {m\ — ml) 2 

r (~9 ->q + q,q + q)=4n f -^ 9 3 q > (10) 

8 m~ 

It is implicitly understood in the subsequent formulae, that all flavors, helicities and C 
related modes for the gluino decays are summed up. Adding up virtual corrections and 
gluon bremsstrahlung, the final result including one-loop SUSY-QCD corrections may be 
written in the form: 

r(g) = r (g){l + ^ [C A F A + C F F F + n f F f + F t + F ren + F dec ]} (If) 

T(g) = T (g) {l + ^ [C A (F A - vr 2 ) + C F (F F + vr 2 ) + n f F f + F t + F ren + F dec ] } (12) 

lr rhe remaining 1/e poles are associated solely with infrared and collinear singularities, while all ultra- 
violet singularities are cancelled. 
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The coefficients Fa etc. are universal functions for both processes, depending only on the 
respective mass ratios 



m§ mj 
r = — | and t = — §■ 
vnA mi 



The difference between the formulae (|TIj) and ( O ) for the two processes is the result of the 
analytic continuation of the universal virtual corrections from the region nig > m g to the 
region m g > rriq, which gives rise to the 7r 2 terms. The coefficients read explicitly: 



r — 1 



Li 2 (l-r) 



r — 1 



Li 2 (-r) + 



5r — 6 



59 



71 



+ 



+ 



3 + r 
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2 



log(r) 
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+ 



r — 1 r — 1 

r 2 r + 1 
r ~ ~2 ~ r - 1 
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log(r) 
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%h2 



Re 



Re 



2 + xi log(l - 1/xi) + x 2 log(l - l/x 2 ) 
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The Fp part represents the supersymmetric version of the abelian vertex and self-energy 
contributions. The Ff and F t contributions are due to quark-squark self-energy corrections 
of the gluino; the gluino-gluon loop is incorporated in Fa- In addition we have introduced 
the renormalization functions 



3Ca — rif — 1 



log 



dec 
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where /Ir denotes the renormalization scale of the coupling a s . The second part Fd ec is added 
to decouple the heavy particles (top quarks, squarks, gluinos) from the running a s . [This 
term may be omitted if the running coupling is adjusted properly]. 

Characteristic features of the corrections are exemplified in Fig. |2]. The three figures 
in the left column display the width q — > gq, the figures in the right column the width 
g — > qq,qq for five light quark flavors. 

The corrections for the squark decays are moderate. In Fig. 0(a) the Born approximation 
(LO) and the next-to-leading order (NLO) SUSY-QCD results are displayed for fiR = rriq. 
The masses are chosen to be irtq = 300 GeV, m- g > 100 GeV, and m t = 180 GeV. The 
running coupling a. s {n R ) is normalized such that a s {M\) = 0.118. As outlined before, the 
running of the coupling is defined in next-to-leading order by the standard QCD spectrum 
with rif = 5 flavors, and the heavy particles are decoupled. The size of the corrections^, 
measured by the ratio Tnlo/^lo m Fig. 0(b), depends on the mass ratio mg/rriq. With 
corrections between +30% and +50%, they appear sizable but still under control. 

The gluino decays are discussed, in a similar way, in Fig. 0(d) and (e). The corrections 
are considerably smaller in this case. The dip in the NLO curve is due to the external 
gluino self-energy, related to the g — > ti threshold in the stop-top loop. The singularities are 
artificial; they are smoothed out by the strong stop-top interactions and the finite lifetimes 
of these particles. Including the widths for t and i in F t removes part of the irregularity (see 
full curves). In view of the small region in parameter space where this effect plays a role, we 
refrain from a more detailed discussion. 

It is evident from the last figures in each column, Fig. Q(c) and (f), in which the variation 
of the widths with the renormalization scale is shown, that the NLO corrections improve 
the stability of the theoretical predictions for the widths significantly. This is the case in 
particular for the gluino decay width, where a wide and shallow extremum is observed for 
I^R/nT'g near unity. 

To shed more light on these somewhat complicated results, it is interesting to consider 
two limiting cases. We have considered the massless limit for the particles in the final states, 
and the threshold effects in the limit of nearly equal masses for squarks and gluinos. For 
the sake of simplicity, we do not decouple the gluinos and squarks from a s and we take for 
illustration also the top mass equal to zero. Setting the scale fiR equal to the mass of the 
decaying particle, we find: 

2 Since Aqcd is not really well-denned in LO, we have identified a s in LO with a s in NLO at the mass 
scale chosen for comparison. 
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massless gluinos/squarks in the final state: 



m 



r (g) <! i + - 
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threshold limit rriq ~ m^: 
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The different size of the corrections for squark and gluino decays is nicely illustrated by 
the limiting case of massless particles in the final state. The driving term for squark decays 
is the large Ca contribution while other color structures play a minor role. For gluino 
decays, however, the different tt 2 terms, induced by the vertex correction for rrig > trig, 
lead to a destructive interference between the reduced (negative) Ca contribution and the 
enhanced (positive) contribution from the other color structures. As a result the SUSY-QCD 
corrections to gluino decays are significantly smaller than the corrections for squark decays. 

In addition to the proof that the effective Yukawa and gauge couplings are identical in the 
limit of exact supersymmetry, we have performed several other consistency checks. We have 
shown, by using a small gluon mass in the DR scheme, that infrared and final-state collinear 
divergences do not give rise to finite differences between the MS and the DR schemes. The 
limiting case rrtg — > discussed in the previous paragraph has been extracted from the general 
formulae and it has been compared with an ab initio rederivation of the squark width (which 
can be calculated much more easily in this limit). Finally, by limiting the calculation to the 
subset of abelian diagrams (i.e. the CV part), we have reproduced the results of Ref. |§ for 
SUSY-QCD corrected squark decays to quarks and photinos after having corrected the sign 
of the vertex correction Eq. (19) in Ref. 0. 
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Appendix 



In this appendix we present the exact supersymmetric one-loop corrections to the qqg and 
qqg couplings, Eq. (||), discussed in Sect. || for Q 2 <C m~ = m 2 q . The corrections are given in 
both the MS and the DR renormalization schemes. 



MS Renormalization Scheme: 

vertex corrections 

a. 



(qqg) 5 V = <J C A 



O's 



(qqg) S v = -r^Cp 
wave-function renormalization 
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wave-function renormalization 



(Z q - l)/2 
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Figure 1: The Feynman diagrams for squark and gluino decays, (a) Born diagram for both 
decay modes; (b) quark, squark, and gluino self-energies; (c) vertex corrections for squark 
decays; (d) gluon radiation for squark decays. The diagrams of type (c) and (d) for gluino 
decays can be generated by rotating the diagrams for squark decays. 
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Figure 2: The SUSY-QCD corrections to squark and gluino decays, (a) The decay rate for 
squarks (rriq = 300 GeV) in LO (dashed curve) and in NLO (solid curve) for ixr = rriq] (b) the 
ratio of the decay rate for squarks in NLO and in LO for ixr = rriq; (c) the scale dependence of 
the decay rate for squarks in LO and in NLO (m g = 200 GeV). (d-f): The same quantities for 
gluino decays with interchanged, but otherwise identical mass parameters. The long-dashed 
curve in (d) and (e) represents NLO results without taking into account non-zero widths for 
stop and top. 
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